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We discuss the renormalisation properties of parity-odd AB = 2 operators with the heavy quark 
treated in the static approximation. Via twisted mass QCD (tmQCD), these operators provide 
the matrix elements relevant for the B° —B° mixing amplitude. The layout of a non-perturbative 
renormalisation programme for the operator basis, using Schrodinger Functional techniques, is 
described. Finally, we report our results for a one-loop perturbative study of various renormalisa- 
tion schemes with Wilson-type lattice regularisations, which allows, in particular, to compute the 
NLO anomalous dimensions of the operators in the SF schemes of interest. 
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1. Introduction 

The oscillations of the system B° — B° are one of the crucial topics in particle physics. Their 
understanding represents a challenging bridge towards the numerical determination of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix and a severe test of the Standard Model. The transition ampli- 
tude responsible for the mixing, 

(BVw+Aa|B°> = \fW B BB, (1-1) 

is mediated by the four-quark operator ^Vv+aa = (py^iby^d) + iby^ys^^y^ysd). It has been 
shown that the renormalisation of such operators is non-trivial in Wilson-like regularisations, re- 
sulting in a mixing with other four-quark operators [[Tj]. Here we propose a strategy to compute 
the matrix element, based on the static approximation of the heavy quark plus the adoption of 
a tmQCD regularisation for the light one. It will be proved that, following these assumptions, 
the mixing under renormalisation is eliminated. Of course, the potential results of the proposed 
approach constitute an intermediate step to the physical solution, as they must be considered in 
view of the calculation of heavy quark subleading corrections and/or interpolations to relativistic 
calculations performed at accessible heavy quark masses [0]. 



2. Operator mapping in tmQCD 

In order to implement our strategy, we start by fixing the notation. The b quark is replaced by 
an infinitely massive quark, described by a pair of static fields (Yh,¥h) propagating forward and 
backward in time, whose dynamics is governed by the Eichten-Hill action ^ (or one of its ALPHA 
variants [Q]), 

S^lYn, Wh) = a 4 £[y h (x)V* ¥h (x) - %{x)V m {x)\ . (2.1) 

X 

On the light quark side, the degrees of freedom are represented by an isospin doublet XjfT = (u,d), 
made of an up and a down quark, and described according to the tmQCD action 1 , 

5 tmQCD = a A £ { ^ ^ + ^ + -^3^] ^ } _ (2 2) 

x 

The equivalence of this regularisation to ordinary QCD, established in [0J, is based on axial trans- 
formations of the quark fields (plus the corresponding spurionic transformations of the mass pa- 
rameters mg and ji( ), which induce a rotation of composite operators between the two theories. In 
particular, for the operator under study one has 

(^vv+A A )f ° = cos(a) (^ v+AA )^ QCD a - isin(a) (^va+av)^'" , (2-3) 

where the terms have to be interpreted as operator insertions in renormalised Green functions in 
the continuum limit, and a mass-independent renormalisation scheme is assumed. Following the 
notation of [B], the twist angle a depends upon the renormalised mass parameters through the 

'We will always work in the so-called twisted basis. For a discussion of the problem in the physical basis, see o]. 
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relation tan(a) = jU^r/w^r, and (2.3) is an identity holding at each value of a. In particular, at 
a = 71 / 2, which is known as the fully twisted case, (2.3) simplifies to 

(<WA)f D = -^VA+Av)?*^ 2 • (2-4) 

In this way, ^Vv+aa i n standard QCD is mapped onto its counterpart &va+av i n tmQCD. Using the 
mass independence of the renormalisation scheme, we will show in the next section that <^Va+av 
renormalises multiplicatively in the static approximation, which represents the main advantage of 
using the above mapping. In this sense, the proposed approach represents an extension to the static 
case of the tmQCD framework used to determine the Bk parameter [Q, 



3. Renormalisation pattern 

We now concentrate on the renormalisation properties of heavy-light four-quark operators, 
with the aim of proving that €?va+av renormalises multiplicatively. Unfortunately, for brevity's 
sake, we skip algebraic details [JTofa . We start by considering generic four-quark operators 

°f ir2 = j Kv^ri^iXv^y^) ± Wiv^O^yi)] , (3.1) 

where Ti 2 represent Dirac matrices. In principle, operators corresponding to different Dirac struc- 
tures could mix among them under renormalisation, thus giving rise to a matrix renormalisation 
pattern; consequently a complete basis of such operators must be considered, such as 



parity — even : Qf 


-o ± 

— ^VV+AA ' 


parity — odd : £}f 


= o% K _ 


fAV 




— ^SS+PP ' 




-o ± 

— ^SP^ 


PS 


Qt 


-o ± 

— I^VV-AA ' 


^3 


= 0^ 


-AV 


ai 


-o ± 

— ^SS-PP ) 




= o±_ 


-PS 



(3.2) 



The renormalisation matrix Z, whose size is in principle 8x8 (mixing between + and — operators 
is trivially excluded), can be constrained through symmetry arguments. Given a symmetry of the 
theory, and the matrix <J> that implements a symmetry transformation at the level of the operator 
basis, it is sufficient to require that Z is invariant under a O-rotation i.e. 

,-1 



Z = <£Z4> 



(3-3) 



The symmetries we use are: 

• Parity. It prevents the mixing among operators with opposite parity. After implementing it, 
the renormalisation matrix Z is reduced to a block-diagonal form, where two 4x4 diagonal 
blocks describe the mixing of the parity-even and parity-odd operators among themselves. 

• Chiral simmetry. It is used a la [|T]] : were chirality respected by the regulator, there would 
be no chance of mixing among different chirality sectors. The mixing due to the Wilson 
chirality breaking in the parity-odd sector can be represented according to the form 







(*n 


<a?± 
°*12 





\ 








ay± 
^22 








0^3 










ay± 
° £ 33 


d¥>± 

"^34 




R 


I 





oy± 
• £ 43 


3?44 ) 



1 + 



/ 





A 31 A 32 



Af 3 Af 4 \ 
A± Af 4 




J 



v4y 



(3.4) 
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where the coefficients J^j are scale dependent, while the Ay's are not. 

Heavy quark spin symmetry and H(3) spatial rotations. We then consider two finite spin 
rotations of the heavy fields, plus two lattice spatial rotations of both heavy and light fields 



heavy quark spin rotations: 
lattice spatial rotations: 



Vh^VhYsYh Wh^WhYiYu 
1 — > 2 ) rotates the 1 axis onto the 2 axis, 

M{ 1 — > 3 ) rotates the 1 axis onto the 3 axis. 



(3.5) 



After a change of basis and some tedious algebra, the parity violating block reduces to 



\£tf - 2=S 4/ / R 



(3^ \ 
0^0 
0^0 

\ 3^ J 



1 + 



fO Af 0\ 
A^ 
Af 
\0 Af oj 



( 4f \ 

=*^3 ~\~ 2JS 4 
\=S 3 — 2J2 4 J 



(3.6) 



• Time reversal. We finally consider a time reversal transformation of the quark fields: 

Vh(x) -> 7o75V/,(* T ), Vh( x ) ^YoYsVhix' 1 ), Wk{x) -» YoYsYki**), k = l,2. (3.7) 

It further constrains the parity-odd block by forcing the residual A,- coefficients in (3.6) to 
vanish. Purely multiplicative renormalisation of J2y A+AV follows therefrom. 

4. Renormalisation in Schrodinger Functional schemes 

We use the Schrodinger Functional (SF) to define a family of finite volume renormalisation 
schemes, in view of a non-perturbative study of the running of the 0va+av operator. Our approach 
closely follows here refs. [f7|], to which the reader is referred for unexplained notation. We first 
introduce bilinear boundary sources at xq = 0,T (being T the time extension of the SF), 



",y x,y 



(4.1) 



where T is a Dirac matrix and the flavour indices si t 2 can assume either relativistic or static values. 
Then, we define a set of SF correlators in order to probe the operators J2f, 



fi lS2 = -jj{KsM^2sAY5}) 



IJIS2 



■^tiKM^M) ■ (4.2) 



The triple fT"^^,^] has to be chosen such that Fj is non-zero. The boundary correlators f\ and 
k\ , which can have either light-light or heavy-light flavour structure, are needed in order to cancel 
the renormalisation of the boundary sources in F l . In practice, we consider ratios of the form 



hf(xo) 



/«[/«]«[*«]!/ 



2- a 



(4.3) 
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and then impose in the chiral limit the renormalisation condition 

5fAf(r/2) = Af(r/2)| w=0 . (4.4) 

Of course, the renormalisation factor depends upon all calculational details, e.g. the light quark 
action (Wilson with(out) a clover term), the static action (Eichten-Hill, or its ALPHA variants), 
the choice of the Dirac structures [Fi,^,^], the value of the 0-angle of the SF and the value 



of the parameter a introduced in ( |4.3[ ). This richness of degrees of freedom can be exploited in 
order to identify some optimal renormalisation schemes, according to the general requirements of 
maximisation of the nonperturbative signal/noise ratio, slowing down of the operator running and 
minimisation of lattice artefacts. 



5. NLO anomalous dimension of £l[ from perturbative matching at one-loop order 

In order to gain information about the running and its lattice artefacts, we have performed 
a one-loop perturbative calculation of the renormalisation factor in some of the SF schemes 
discussed above. Such a calculation allows us to determine the NLO anomalous dimension of £2^ 
via a perturbative matching to some reference scheme in which the NLO anomalous dimension 
is already known. The matching procedure has been illustrated and applied several times in the 
literature [ |T"i| , [7[ [TJ|], and it will not be reviewed here. The reference scheme was chosen to be 
DRED, where the NLO anomalous dimension of and its perturbative matching to the so-called 



to-scheme have been computed in [jl3j]. The perturbative expansion of reads 



^ + (gl,a/L) = l + Zgf^ k) = l+g 2 Q [tf + tf kg) +0(l)\ +0(gt) (5.1) 

k=l 

where = —l/(2n 2 ) is the universal anomalous dimension of the operator £}f, and is the 
one-loop scheme-dependent finite part, peculiar to the SF and the defining choices listed at the end 
of the previous section. The running of the operator is described by the step scaling function (ssf), 



cij («) = hm 



= 1 + E«8V*> (5-2) 

g 2 (L)=u k=l 



As an example of the running, the ssf of at NLO and Nf = is reported as a function of the 
renormalised coupling on the left side of Figure 1. The plot refers to the choice fTi,^,^] = 
[75) 75) 7s]- The straight line represents the universal LO running, and the bands describe the depen- 
dence of the NLO anomalous dimension upon the choice of a, when the latter ranges in the interval 
[0, 1 /2\. On the right side of Figure 1 we report a comparison of the lattice artefacts 8^{a/L) on 
the ssf, defined as in [ 1 1], between the static-light case and the light-light one (data from [§]). The 



comparison refers to the schemes where [ri,r2,r3] = [75,75,75], 6 = 0.5 and a = 0. The light 
quarks are discretised according to the unimproved (W) or the c sw -improved (SW) Wilson action, 
while the static quarks are discretised according to the Eichten-Hill (EH) action. Although the 
static-light schemes cannot be directly compared to the relativistic ones (where the normalisation 
of the four-quark correlator is always performed using only the relativistic correlators), the plot 
shows that the introduction of static quarks does not imply a significant increment of the lattice 
artefacts in perturbation theory. 
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g 2 a/L 

Figure 1: On the left side the step scaling function of at NLO and Nf — is reported vs. the renormalised 
coupling in the SF scheme. On the right side we compare the lattice artefacts of the step scaling function 
between the full relativistic case and the static-light case. Both plots are preliminary. 
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